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Keeping in mind the possibility of large #13, which will be soon ex- 
plored by reactor and accelerator experiments, I formulate a perturbation 
theory of neutrino oscillation under the ansatz s 2 3 ~ Am| 1 /Am§ 1 ~ 0.03, 
which is comparable to the Chooz limit. Under the framework, I derive the 
perturbative formula of the v e appearance probability valid to order e 2 in 
which effects of arbitrary matter density profile is taken into account. I use 
the formula to analyze problem of possible obstruction to detecting lepton 
CP violation by effects of asymmetry in matter density profile. Though 
the asymmetry could be large for neutrino trajectories which traverse both 
continental and sea crust, its effect on obscuring CP violation measurement 
is found to be quite small. 

PACS numbers: 14.60.Pq,14.60.Lm,91.35.-x 

1. Introduction 

My presentation in Ustron'09 was done under the title "Long-Baseline 
(LBL) Neutrinos; Looking Forward to the Future" . It included a review of 
the ideas for exploration of the unknowns in the 1-3 sector of the MNS ma- 
trix pQ , CP violation due to the lepton analogue of the Kobayashi-Maskawa 
phase \'2\ 5 and the neutrino mass hierarchy. If I restrict myself into per- 
spective in North-East Asia, they include an upgrade of J-PARC beam with 
megaton-scale Hyper-Kamiokande (T2K II) as described in [3], a 100 kt scale 
liquid Ar detector in Okinoshima [4J, and the Tokai-to-Kamioka-Korea (for 
short T2KK) setting [51 [6]. The overview of the latter is given in [7]. In 
particular, I emphasized the Kamioka-Korea identical two-detector setting 
as a robust way of measuring the CP violating phase and determining the 
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mass hierarchy. It combines the idea of low energy superbeam as the clean- 
est way for detecting lepton CP violation [8] and the powerfulness of the 
two-detector method [9]. But, since the contents of this part are described 
in the previous reports |10t [TT] I confine myself into the last part of my 
presentation in Ustron'09, the large-#i3 perturbation theory in this written 
version. 

2. Motivation and Use 

The motivation for formulating the large-#i3 perturbation theory is very 
simple; #13 can be as large as the Chooz limit [12]. If it is the case S13 ~ 0.17. 
I emphasize that this possibility is to be tested very soon by the accelerator 
[3j [13] and the reactor #13 experiments [2] some of which will start in 
2009. Then, the e perturbation theory (in a terminology defined in [15j ) 
formulated under the ansatz S13 ~ e = Am^i/Am^ — 0.03, which provides 
the simplest way of deriving the widely used Cervera et al. formula [16] of 
the oscillation probability, would not serve as the best approximation. 

Then, what is the use of the large-#i3 perturbation theory? I confine 
in this paper the robustness issue in uncovering lepton CP violation. As is 
well recognized the matter effect produces a fake CP violation. The presence 
of the matter effect is inevitable for settings which also have sensitivity to 
the mass hierarchy. Early references of this topics include |1 7|. 118] . The 
problem has been discussed in a number of authors which produced too 
many references to quote here. However, there exists a point which does not 
appear to be given full attention in the literature, the problem of possible 
asymmetry in matter density profile in the earth. See, however, [19] for 
discussion of this problem. It is known that asymmetric baseline produces 
CP violating sin J terms in P(un — > Vu) [2D], and cos 5 terms in T violating 
observable APt [21] which invalidates the neat property of T violation, 
APt = for vanishing 5 [22] . Apparently, the effects of asymmetric matter 
density profile is most prominent for large #13. 

A crucial question would be: Are there any situation in which sizable 
asymmetry in matter density profile shows up? The answer is indeed yes. 
Suppose a baseline of order ~ 1000 km and a neutrino beam launched at 
a place on a continent is received by a detector which is placed in an is- 
land. If the travel distances of the beam in the continent and in the sea 
are comparable, one can expect asymmetry in the density profile because 
the matter density in the earth crust is believed to be smaller under the 
sea. Contrary to the mantle density, the crust density is not severely re- 
stricted by the earth mass, and it would not be easy to measure it directly 
either. Therefore, it is important to investigate to what extent CP violation 
discovery might be obscured by asymmetry in matter density profile. 
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3. Large-#i3 perturbation theory 

The neutrino evolution equation can be written in flavor b 



H a pv/3 (a, p = e,fi, t) 
tonian is given by 



dx 
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In the standard three-flavor neutrino scheme, Hamil- 
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where Am^ = rnj — mf, and a(x) = 2\^2GpN e (x)E is the coefficient which 
is related to the index of refraction of neutrinos in medium of electron 
number density N e (x), where Gf is the Fermi constant and E is the neutrino 
energy. U = U23U13U12 is the MNS matrix [1] in the lepton sector. 

It is straightforward to formulate the perturbative framework of neutrino 
oscillations. We refer [15] for notations. We use the tilde-basis v 
with the tilde-Hamiltonian H = U23HU23, which is decomposed as H 
Hq + Hi. Then, the S matrix can be written as 



S{L) = U 23 e- iHox n(x)U. 



; 23- 



(2) 



Q(x) can be expanded with use of Hi = e tH ° x Hie tH ° x as 

Q{x) = 1 + (-i) f X dx'Hi(x') + {-if f X dx'Hi(x') [* dx"H x {x") + 0(e 3 X3) 
J Jo Jo 

where the "space-ordered" form in ([2D is essential because of the highly 
nontrivial spatial dependence in Hi. 

We use the method of Fourier decomposition to incorporate the effect 
of matter density variation in the earth [23] with the dimensionless variable 
v a = A ^ ■ It can be expanded into a Fourier series as 



n=l 



r A e ~ip n x _|_ ^ r Ay e ip n x 



(4) 



where p n = ^n. It should be noticed that if ta{L — x) = ta(x), namely if 
the baseline is symmetric, = (r„ )*. Therefore, the imaginary part of 
represents the effect of asymmetric matter density profile. 

We assume large #13 comparable to the Chooz limit, S13 ~ y/e, where 
e = Amli/Amfi — 0.03. Other small parameters are present depending 
upon experimental setting, in combination with matter density, neutrino 
energy, and baseline. As a model superbeam experiment we consider a 
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baseline with distance L 
Noticing that -^V = 0.084 

Ami 



1000 km and neutrino energy E 

a 2 \ - 1 

2.5xlO~ 3 eV / 



(idlv) 



2.8g/cm" 



2 GeV. 
, the set- 



ting leads to a/Am^ = 0.17 ~ y/l. (For T2K II setting, a/Am^ ~ e 
may be more appropriate.) Then, we formulate the large-#i3 perturbation 
theory by taking the following expansion parameters e = Am|i/ Am^ and 
s i3 — a/Amjj ~ y/e. The unperturbed part of the tilde-basis Hamiltonian 
is given by Hq = Adiag(0, 0, 1), while the perturbed part is 
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where A = Am 2 1 /2£'. The first, second, third, and the fourth terms in ((SJ) 
are of order ez, e 1 , ea, and e 2 , respectively. 

The S matrix can be computed with use of Then, the appearance 
oscillation probability P{yn —> v e) 



with a simplified notation A31 = Am^L/AE by 



|5 eM | is given to second order in e 
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4. T, CP and CPT violation observable 

It is interesting to compute the expression of T, CP and CPT violation 
observable to know the difference in their dependence on 5 and the matter 
density profile. The T, CP, and CPT conjugate probabilities can be obtained 
by the appropriate replacements: 
P(ve -> = P -> v e ; -5, t$, Re(r^), -Im(r^)) , 

Pipp -^u e )=P fa -> v e ; -5, -r$, -Re(r^), -Im(r^)) , 

P{v e -> ^) =-P(^ -> ^ e ;<5,-r^,-Re(r^),Im(r^)). 

With © it is it is easy to calculate them and the results read: 

T violation : Apr = P{ya — ► v e) — P{y e — > v^) 
= — 16 J r sin 5eAsi sin 2 A31 

- 16 J r sin SeA 3i (sin 2 A 31 - - A 31 sin 2A 31 j r£ 



- 32 J r sin 5eA 31 sin 2 A31 £ ^2 Re (rf) 

n=l A 31 V n7r i 

00 A 2 / \ 

- 32 J r cos 5eA 2 x sin 2 A31 ]T f 2 31 7 . 21 Im (r£) 

00 

+ 32 S 2 3S 2 3 A| 1 sin 2 A 31 r ^ £ ™ Im (tf) . (7) 

n=l ^31 ~~ l n7r J 

CP violation : AP^p = P^ — ► i/ e ) — P(^t — > i? e ) 
= — 16 J r sin 5e A 3 i sin 2 A31 

+ 16s^ 3 s 2 3 (sin 2 A31 - ^A 31 sin2A 31 ^ r^ 1 

OO A 2 

+ 3244 sin 2 A 31 £ A 2 1^2 Re 

n=l 31 V / 

+ 16 J r cos 5eA 31 ( A31 sin 2 A 31 + - sin 2A 31 - A31J r£ 



+ 16 J r cos 5eA 31 sin 2A 31 £ ^^Re (r£) 

n=l 31 V n7F J 
°° A 2 / \ 

- 32 J r cos JeA 2 ! sin 2 A31 £ r 2 31 , , 21 Im (r£) . 

CPT violation : APqpt = ^(^u ~~ ► v e) ~ Pipe ~~ * ^u) 
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= I6S23S13 ^sin 2 A31 - ^A 3 i sin2A 3 i^ rfi 

OO A 2 

+ 32^4 sin 2 A 31 E A , ^ H 1 *) 



1 "31 



+ 16J r eA 3 i 



A31 sin A31 sin (5 + A31) + sin A31 cos (5 + A31) - A31 cos 5 



'0 



00 ^2 

+ 32 J r eA 3 i sin A 31 cos (5 + A 31 ) £ ^—^ Re (r^) 

n=l 31 i n7r J 

00 

+ 3244AI, sin 2 A 31 r ^ £ "1^ *" «) . (9) 

n=l ^31 ~~ V n7r J 

In (JJJ), (jHJ), and ([9]), first a few terms are of order e 3//2 , and the rests are of 
order e 2 . The effects of asymmetric profile are always contained in the latter. 
As we noted earlier the charming property of T violation observable holds 
without asymmetry in the matter profile, Im(r„ ) = 0. That is, if 5 vanishes 
then AP? = 0; The matter effect cannot produce a fake T violation. The 
asymmetric baseline destroys the neat property. 

We examine the effect of asymmetry in the matter density profile by 
taking an explicit model of the profile: 

p = p + 5p (0<x< L/2), 

p = p -5p (L/2<x<L), (10) 

which leads to lm(r^) /r^ = (2/im)5p/ po (odd n), and Re(r^) = 0. r A = 
ao/Am 2 ! where ao is obtained by using the matter density po in the defini- 
tion of a. For concreteness we take po = 2^/cm 3 and bp = 0.8g/cm 3 . We 
use the toy model of matter density profile to compute T, CP, and CPT 
violating observable APt, APqp, and APqpt defined respectively in (JT]), 
dU) , and © . We do this under the approximation of keeping only the low- 
est mode Im(r^) (n = 1), which appears to give a good approximation. I 
take sin 2 2# 13 = 0.1 and 5 = 3tt/4 in this calculation, but it appears that 
qualitative features are rather insensitive to 5. 

The results of the calculation are presented in Fig. [TJ One can see 
immediately that, in spite of a rather large asymmetry in density profile 
p0|) the effect of the asymmetric baseline (shown by the red lines) is small. 
In fact, it is negligibly small in CP violation observable. It appears that 
this feature prevails for other choices of CP violating phase 5. I have also 
checked that the effect of the asymmetry becomes even smaller for shorter 
baseline and/or smaller #13. I believe that this settles the issue of possible 
contamination effects by the asymmetric baseline to detection of genuine CP 
violation effect; A good news for medium baseline (~ 1000 km) experiments. 
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Fig. 1. T, CP, and CPT violating observable AP T , AP C p, and AP C pt defined 
respectively in ([7]), ([5]), and © are plotted as a function of neutrino energy (in 
GeV) in the upper, middle, and the lower panels, respectively. In each panel the 
blue and the red lines indicate, respectively, the value of AP and the contribution 
to AP from Imfjj 4 ) terms. The green lines in the upper and the middle panels 
indicate the contribution of vacuum effect only. The green line in the bottom panel 
represents the contribution from average matter density. 
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Some remarks are in order: 



• The size of the matter effect (difference between blue and green lines) 
gives a dominant effect in APqp, overturning the negative sign of 
the vacuum contribution in this particular case. But, it plays only a 
minor role in APt- CPT-violation observable would be most powerful 
to resolve the mass hierarchy [26J because it gives the largest AP. 

• In CP (also T) violation observable, the energy dependence of the 
average matter density term (blue minus (green + red)) are rather 
similar to the vacuum term. Then, there could be severe confusion 
between CP violation caused by phase and uncertainty in the average 
density of matter. In particular, it gives rise to a serious confusion 
at 5 = because of its cos 5 dependence. Careful spectrum analysis 
would be required to resolve the confusion between the matter-CP and 
phase-CP effects. 

• For smaller S13 ~ e, which may be relevant for neutrino factory set- 
ting, one can show by using the e perturbation theory [15] with small 
matter density variation of ~ e that the terms sensitive to den- 
sity variation are at most of order e 3 . It confirms qualitatively the 
conclusion reached in \23\. 
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